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The Pioneer anomaly is a small sunward anomalous acceleration found in the trajectory analysis 
of the Pioneer 10 and 11 spacecraft. As part of the investigation of the effect, analysis of recently 
recovered Doppler data for both spacecraft has been completed. The presence of a small anomalous 
acceleration is confirmed using data spans more than twice as long as those that were previously 
analyzed. We examine the constancy and direction of the Pioneer anomaly, and conclude that: i) the 
data favor a temporally decaying anomalous acceleration (~2x 10~ n m/s 2 /yr) with an over 10% 
improvement in the residuals compared to a constant acceleration model; ii) although the direction 
of the acceleration remains imprecisely determined, we find no support in favor of a Sun-pointing 
direction over the Earth-pointing or along the spin-axis directions, and iii) support for an early 
"onset" of the acceleration remains weak in the pre-Saturn Pioneer 11 tracking data. We present 
these new findings and discuss their implications for the nature of the Pioneer anomaly. 

PACS numbers: 04.80.-y, 95.10.Eg, 95.55.Pe 
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I. INTRODUCTION 

Analysis of the navigational tracking data received 
from the Pioneer 10 and 11 spacecraft at large helio- 
centric distances of ^20-70 AU indicated the presence 
of a small anomalous Doppler frequency drift in their 
radio- metric observables [l|, [2| ■ Ultimately, this drift was 
interpreted as an anomalous constant acceleration act- 
ing on both of these spacecraft in the sunward direction, 
with a magnitude of a P = (8.74 ±1.33) x 10" 10 m/s 2 Q. 
There were also earlier reports that, in the case of Pioneer 
11, the anomaly may have begun with a relatively sud- 
den "onset" [3| shortly after the spacecraft's encounter 
with Saturn. This acceleration of unknown origin is to- 
day known as the Pioneer anomaly (for a review, see |4|). 
In this paper, we present an analysis of newly recovered 
Doppler data for both Pioneer 10 and 11, and probe the 
question of the constancy of the anomalous acceleration, 
as well as place constraints on its direction. 



II. DATA SETS AND ANALYSIS 

Following the 2002 study 3], an effort was initiated 
to collect all available Doppler data for both spacecraft 
d, Q . Using standard data conditioning techniques [3| , 
several decades-old archival data tracking files (see dis- 
cussion in |5j) were converted to modern formats, pro- 
cessed, edited and filtered for corrupt data. Uplink fre- 
quency records at NASA's Deep Space Network (DSN) 
were reconstructed from redundant information. We 
summarize the newly recovered Pioneer 10 & 11 Doppler 
data sets in Table HI including the range of dates covered 
and number of two-way and three-way [3j, |4| coherent 



Doppler data points available. The largest portion of the 
data lies in "deep space" ("DS") beyond planetary en- 
counters (heliocentric distances of 18-80 AU and 9-32 
AU for Pioneer 10 & 11 respectively). The data also in- 
clude planetary encounters with Jupiter ( "J" ) and Saturn 
("S" , for Pioneer 11 only). Data for other time ranges did 
not survive to the present. Here, we focus on the deep 
space data, and the portion of Pioneer 11 Saturn data 
which is primarily under the gravitational influence of 
the Sun, which we designate Saturn "approach" ( "SA" ) . 
These latter data are outside of Saturn's Roche (or Hill) 
radius so that Saturn and its moons can be treated only 
as a perturbing influence. 

Compared to the 2002 data, the data arc for Pioneer 10 
has doubled in length, from 11.5 years to 23.1 years, and 
the total number of data points increased from 20,055 
to 41,054, including Jupiter encounter. For Pioneer 11, 
the length of the contiguous data arc increased from 3.75 
to 10.75 years, and the number of data points increased 
from 10,616 to 81,537, including encounters. In terms of 
heliocentric distance, the Pioneer 10 & 11 arcs overlap 
from 18-32 AU. 

During the mission, both spacecraft performed numer- 
ous maneuvers. As the Pioneers were spin-stabilized, pe- 
riodic attitude maneuvers were required to realign the 
spin and antenna axis to the spacecraft-Earth direction 
to within 1.5°. While attitude maneuvers used a fore-aft 
pair of thrusters in tandem, small residual impulses along 
the spin direction may occur for [4j. Information about 
maneuvers is available from recovered mission records 
and from the spacecraft telemetry [5j, and is summa- 
rized in Table |TJ The data given are lower limits as it 
is possible that a few maneuvers are not reflected in the 
available telemetry. 



TABLE I: Pioneer Doppler Data Sets 



s/c a 


Data 6 


Date 


Range 


Points 


Man. c 


P10 


J 


1973-10-15 


- 1973-12-27 


5806 


>6 




DS 


1979-02-14 


- 2002-03-03 


35248 


>83 




2002DS 


1987-01-03 


- 1998-07-22 


20055 


28 


Pll 


J 


1974-04-18 


- 1974-12-26 


7467 


>16 




S 


1977-11-01 


- 1979-09-18 


9017 


>35 




SA 


1977-11-01 


- 1979-06-29 


4282 


— 




DS 


1980-01-12 


- 1990-10-01 


65053 


>92 




2002DS 


1987-01-05 


- 1990-10-01 


10616 


22 



a (P10), (Pll) denote Pioneer 10 and Pioneer 11 respectively. 
'(J) Jupiter; (S) Saturn; and (SA, subset of S) Saturn approach; 
(DS) deep space; (2002DS) is the data used in the 2002 study 0. 
c Number of attitude maneuvers performed. 



Our analysis was carried out using the Jet Propul- 
sion Laboratory's (JPL) Orbit Determination Program 
(ODP) that was used for earlier work on the anomaly 
[3,|j]. As before, central to the analysis was establishing a 
model orbit for the spacecraft that takes into account all 
known forces, gravitational and non-gravitational, while 
numerically integrating the appropriate equations of mo- 
tion. The model included the effects of planetary per- 
turbations, solar radiation pressure, propulsive maneu- 
vers, general relativity, and bias and drift in the Doppler 
observable. Planetary coordinates and the solar sys- 
tem masses were obtained using JPL's Export Planetary 
Ephemeris DE421. The model also included the precise 
positions of Earth-based stations of the DSN as well as 
radio propagation effects (see details in |J|). 

The study of the deep space data arcs (DS, see TableHJ 
from the extended Doppler data set was expected to im- 
prove our understanding of two key characteristics of the 
Pioneer anomaly: its direction and its temporal behav- 
ior. Additionally, it was expected that some of the early 
data, in particular the Pioneer 11 Saturn approach data 
(SA), would help us confirm whether or not the anoma- 
lous behavior began with a relatively sudden "onset" [3| . 

We considered three models for the anomalous accel- 
eration — constant, linear and exponential — all applied 
along the nominal Earth-spacecraft line. The constant 
model has one parameter, ap, representing a constant 
modeling error. The linear model, 



a P (t) = a P (t ) + (t- t )a P 



(1) 



contains a jerk term, ap. The exponential model, 

a P (t) = ap(to)e-^-^ ln2 (2) 

decays with half life /3 _1 . This last model is physi- 
cally motivated by a potential relation to the on-board 
power generators, which radioactively decay. The epoch 
is t ^January 1, 1972. 

Separately, the anomalous acceleration was estimated 
using a batched stochastic model 0, El- This method 
produces a smoothed acceleration for each batch 3] ■ As 
the model with the most estimated parameters, it is likely 



to produce the best possible fit, but since it is purely 
phenomenological, it provides the least physical insight 
into the anomaly. 

To consider the direction of the anomaly, we separately 
modeled the acceleration as a constant vector in four 
principal directions [3, Q: that of i) the Earth, ii) the 
Sun, iii) the spin axis, and iv) the spacecraft velocity 
vector. Of these, the spin and Earth-spacecraft axes are 
effectively degenerate as the spacecraft were maintain- 
ing an Earth orientation for continuous radio communi- 
cation. We estimated an acceleration vector that was 
constant in a reference frame with its z-axis aligned with 
the spacecraft-Sun line, and an acceleration vector that 
was constant in a reference frame with its z-axis aligned 
along the spacecraft-Earth line. The x-axis in both cases 
lay in the plane of the ecliptic, and the y-axis completed 
a right-handed triad. ODP solves for the components of 
the acceleration vector independently. 

For all acceleration models, the initial state vector and 
the velocity impulses for each attitude maneuver were 
also estimated using a least squares fit. 

The estimate and formal errors for the model param- 
eters are computed by the least squares fit assuming the 
Doppler measurement errors are uncorrelated. In real- 
ity, the residual still shows significant structure, perhaps 
due to mismodeling (e.g., models of solar plasma, the at- 
mosphere, etc.). The presence of such autocorrelation in 
the residuals is the reason why it is common for "real- 
istic" errors to be larger than formal errors by an order 
of magnitude or more. In this Letter we report formal 
errors, with the understanding that these values do not 
necessarily represent the actual uncertainty of model pa- 
rameters, but can be used as an indication of the overall 
goodness of the fit and thus, the quality of the model. 

Acceleration models are compared computing ov, the 
root mean square (RMS) of the difference between mod- 
eled and observed Doppler frequencies. A typical "good" 
fit should yield RMS residuals of 10 mHz or less, often 
under 5 mHz. An RMS residual above ~10 mHz usually 
indicates that the model does not adequately describe 
the physics of the spacecraft's motion, the dynamics of 
the solar system, or the propagation of the radio signal. 
We find that the Pioneer 10 data set is noisier than that 
for Pioneer 11, with RMS residuals differing by a factor 
of two. 

We evaluated the models described in the previous sec- 
tion using the ODP. The deep space data sets of both 
the Pioneer 10 and 11 spacecraft 1 were evaluated us- 
ing the one-dimensional (1-D) constant, linear, and ex- 
ponential models, and the results are shown in Table 
III1 The fit quality of the different temporal models is 
nearly the same. The variable models are consistent with 
a gradually decreasing acceleration, either at a rate of 



1 For technical reasons, Pioneer 10 data up to 1998 and Pioneer 
11 data for 1983-1990 were used. 



TABLE II: Deep Space Acceleration Parameters 



s/c a 


Model 


mHz 


a P b 
lO" 10 m/s 2 


Additional 
parameter 6 ' 


P10 


Constant 

Linear 

Exponential 


4.98 
4.60 
4.58 


8.17(2) 
11.06(8) 
12.22(16) 


a P = -0.17(1) 

/T 1 = 28.8(0.7) 


Pll 


Constant 

Linear 

Exponential 


3.67 
2.09 
2.06 


9.15(7) 
11.65(42) 
13.79(62) 


a P = -0.18(3) 
/T 1 = 24.6(2.4) 



"See Table U for designations. 

b l-a formal error quoted in final digit(s). 

c ap in 10 -10 m/s 2 /yr, [3~ 1 in yr. 



TABLE III: Vector Acceleration Parameters 



S/C a 


Data 


Center 


mHz 


a P J 

10 


b 

ap,x 
" 10 m/s 2 


b 

ap,y 


P10 


DS 


Earth 
Sun 


4.98 
4.98 


8.17(2) 
8.17(2) 


0.42(3) 
0.42(3) 


0.08(1) 
0.08(1) 


Pll 


DS 


Earth 
Sun 


2.13 
2.21 


8.62(6) 
8.64(6) 


-0.73(3) 
-0.76(3) 


0.97(1) 
0.97(1) 


Pll 


SA 


Earth/ 
Sun c 


3.61 


4.48 
±0.50 


-1.38 
±11.33 


-0.27 

±1.7 



"See Table U for designations. 

6 l-cr error quoted in final digit, except as noted. 

c Results for Earth-/Sun-pointing were similar. 



~1.7x 1CP 11 m/s 2 /yr, or a half life of ~27 yr averaged 
for both spacecraft. 

The results of the exponential model are shown in Fig- 
ure [TJ The figure also shows the stochastic estimate, il- 
lustrating the similarity between the parameterized and 
stochastic models. The RMS residuals of the stochastic 
fit show further improvement compared to other models: 



opio = 4.40 mHz, 



o-pu = 2.02 mHz. 



(3) 



The results of the three-dimensional (3-D) models are 
shown in Table IIII1 For the DS data sets, acceleration 
vectors that are constant in a solar system barycentric 
reference frame or constant relative to the spacecraft- 
Earth line cannot be distinguished. The estimated vector 
direction is within 6° of the spacecraft spin axis. 

We also investigated the possibility of an "onset" of 
the anomaly using Pioneer 11 Saturn Approach ("SA") 
data from 1977-1979, the results of which are also shown 
in Table IIIII While we attempted to fit a simple 1-D 
Earth pointing constant acceleration, results were poor 2 . 
The vector model has very large errors in the x direction, 



2 The validity of the SA results was questioned because of the large 
uncertainty in the acceleration estimates, their sensitivity to so- 
lar radiation pressure parameters, the structure of the residual 
and presence of a 21 m/s maneuver. In the vicinity of this ma- 
neuver, the residual display shows residuals of 20 mHz, five times 
greater than the typical residuals outside this region. 



10 



n +16 

x 

A -16 



-8 r 



- 1 1 1 - 



-, — i — i — i — i — |— 
Exponential 



H 1 1 H 



H 1 1 H 



H 1 1 H 



Stochastic 




■#-^f-^-#¥#ii- 



10 15 20 

Years after Launch 



25 



FIG. 1: Top panel: Estimates of the anomalous acceleration 
of Pioneer 10 (dashed line) and Pioneer 11 (solid line) using 
an exponential model. Second panel: Stochastic acceleration 
estimates for Pioneer 10 (open circles) and Pioneer 11 (filled 
circles), shown as step functions. Bottom two panels: Doppler 
residuals of the stochastic acceleration model. Note the dif- 
ference in vertical scale for Pioneer 10 vs. Pioneer 11. 



which relate to a larger solar radiation pressure effect at 
these distances. The acceleration magnitude for Saturn 
approach is a P n = (4.58 ± 11.80) x 10" 10 m/s 2 . This 
result is consistent with the acceleration estimates of the 
DS phase, and thus, we cannot conclude definitively that 
an "onset" exists or not (see [3, |4| for discussion) . 



III. DISCUSSION 

Temporal behavior: We can unambiguously confirm the 
presence of an anomalous acceleration the recently recov- 
ered Pioneer 10 & 11 spacecraft tracking data, with con- 
sistent magnitudes between the two spacecraft, and also 
consistent with previous results [l|, y, la, 0] • 

The constant acceleration models, both 1-D and 3-D, 
are essentially identical for Pioneer 10. However, for Pio- 
neer 1 1 , the RMS residuals improve when considering an 
unknown constant force, perpendicular to the spacecraft- 
Earth direction. This may be related to the anomalous 
spin- up of Pioneer 11 5]. 

Our estimates of a jerk term ((Tl) are consistent with 
earlier studies J4J. Markwardt [a] obtained an im- 
proved fit of Pioneer 10 data when estimating a jerk 
of apio = -0.18 x lO" 10 m/s 2 /yr; also Toth Q ob- 
tained a P io = (-0.21 ± 0.04) x 10~ 10 m/s 2 /yr, a P11 = 
(— 0.34±0.12) m/s 2 /yr for Pioneer 10 & 11, respectively. 

The rationale for an exponential model ^/ is based on 
the possibility that the acceleration may be due to ther- 
mal recoil forces generated on-board. Due to degradation 
of the RTG thermocouples and changes in the thermal 
louver system 4] , the resulting thermal recoil force could 
have a half-life significantly shorter than the 87.74 year 



half-life of the 238 Pu fuel [8], with 27 years being in the 
acceptable range. 

The gradually decreasing linear and exponential decay 
models yield marginally improved fits when compared to 
the constant acceleration model, as does the stochastic 
model. The presence of maneuvers confounds our ability 
to detect such terms unambiguously. The addition of 
earlier data arcs, with greater occurrences of maneuvers 
did not help as much as desired. 

Our measure of goodness-of-fit, <j v , may allow com- 
parison of competing models using the standard F-ratio 
statistic (i.e. ratio-of- variances), even if <j v itself is not 
a standard statistic. Preliminary results indicate that 
the time variable models (linear, exponential, stochastic) 
are indeed better than the simple constant models, for 
both spacecraft trajectories. Future work will attempt 
to quantify this improvement more rigorously. 

Direction: For continuous communication, it was nec- 
essary to orient the spacecraft so as to keep the Earth 
with the 3° of their antenna beamwidth. The Sun-probe- 
Earth (SPE) angle remained small and varied only from 
2.6° (1980) to 0.7° (2001) for Pioneer 10 and from 6.0° 
(1980) to 1.7° (1990) for Pioneer 11. The solar plasma 
noise present in the data rendered the effort to distin- 
guish these nearly coincident directions fruitless. There 
was also a possibility that Pioneer 11 data from prior to 
Saturn encounter, when the SPE angle was larger, would 
allow us to distinguish between the Earth (and spin axis) 
vs. Sun directions. Unfortunately, these hopes were, too, 
in vain because of solar plasma noise, a malfunctioning 
thruster, and frequent maneuvers. Although the Earth 
direction is marginally preferred by the solution (see Ta- 
ble I!!!]), the Sun, the Earth, and the spin axis directions 
cannot be distinguished. 

We can exclude an anomaly directed along the space- 
craft velocity vector. In 1980, the angle between the best- 
fit acceleration vector and the spacecraft velocity vector 
is 8.5° (Pioneer 10) and 31.8° (Pioneer 11). This was suf- 
ficient to show that the anomaly is not directed along the 
velocity vector using Toth's orbit determination program 
Q with the Pioneer 11 DS data arc. 

Onset: The Doppler data obtained using Pioneer 11 
prior to its Saturn encounter are consistent with a pos- 
sible onset, but the uncertainty remains large. Some, 



or all, of this onset may be due to mismodeling the ef- 
fects of solar pressure. However, the relative shortness 
of these data arcs and the large number of maneuvers 
performed during this period (including a significant tra- 
jectory correction maneuver) make it difficult to reach a 
robust conclusion. 

Origin: The most likely cause of the Pioneer anomaly 
is the anisotropic emission of on-board heat. This fact 
was recognized early on [3j, leading to a detailed ther- 
mal analysis of the Pioneer 10/11 spacecraft. In a par- 
allel effort, using recovered project documentation and 
telemetry records, a highly detailed finite-element ther- 
mal model of the two spacecraft was constructed and 
used to estimate the recoil force due to anisotropically 
radiated on-board generated heat at various heliocentric 
distances 0, Q . A conclusive result can only be reached 
by incorporating the thermal recoil force, computed as a 
function of time, into the standard set of spacecraft force 
models that are used for Doppler analysis 4, 9j. Such an 
analysis was initiated once the extended Pioneer 10 and 
11 Doppler data sets became available. The main ques- 
tion is whether or not a statistically significant anoma- 
lous acceleration signal still remains in the residuals after 
the thermal recoil force has been properly accounted for. 
Results of this meticulous study will be published soon. 
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